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Abstract 

The dynamical equations which are basic for the description of the dynamics of quan- 
tum felds in arbitrary space-time geometries, can be derived from the requirements of 
a unique deterministic evolution of the quantum fields, the superposition principle, a 
finite propagation speed, and probability conservation. We suggest and describe obser- 
vations and experiments which are able to test the unique deterministic evolution and 
analyze given experimental data from which restrictions of anomalous terms violating 
this basic principle can be concluded. One important point is, that such anomalous 
terms are predicted from loop gravity as well as from string theories. Most accurate 
data can be obtained from future astrophysical observations. Also, laboratory tests like 
spectroscopy give constraints on the anomalous terms. 

1 Introduction 

Experimental Quantum Gravity, the experimental search for deviations from Einstein's Gen- 
eral Relativity, which includes also Special Relativity, has been an very active area in physics 
since a few years. All unifying theories or quantum gravity theories predict small modifica- 
tions from General Relativity, for example deviations from the 1/r-potential, and violation 
of the equivalence principle, see e.g. violation of Lorentz-invariance and violations of 
the universality of the gravitational red shift. Violations of these basic principles underlying 
General Relativity go together with a modification of the equations of motion for test matter 
in the gravitational field. 

Questions about the structure of the field equations for quantum objects came up very 
recently in the context of quantum gravity: From loop gravity Q as well as from string theory 
H there are predictions about quantum gravity induced modifications of the field equations 
governing the motion of spin-i-particles. Up to now the predictions consist in directional 
derivatives and higher-order spatial derivatives added to the usual Dirac equation. However, 
more general modifications can also be expected. For example, the prediction of second 
order spatial derivatives, which is worked out in a distinguished frame of reference, makes it 
reasonable to expect also second order time derivatives. Therefore, it is important to study 
experimental consequences of general non-standard modifications of the Dirac equation. In 
order to have a guiding principle at hand which tells us something about the meaning and 
the physical consequences of certain modifications of the usual Dirac equation, we will study 
these modifications in the frame of a constructive axiomatic derivation of the Dirac equation. 

In that approach it is possible to derive the Dirac equation from a very few fundamental 
operational principles which can be tested directly in experiments. The derivation of the 
Dirac equation can be divided into two parts: The first part consists in the derivation of 
a system of linear hyperbolic field equations, which we call a generalized Dirac equation, 
from four fundamental principles which the quantum field has to obey, namely (i) the unique 
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deterministic evolution, (ii) the superposition principle, (iii) finite propagation speed, and 
(iv) conservation of probability: 

unique deterministic evolution "] 

superposition principle J generalized Dirac-equation 

finite propagation speed | \0 = i"/ a (x)d a ip(x) — M(x)tp(x) 

conservation law J 

Here 7 (a = 0, . . . ,3) are some matrices which are not assumed to fulfill a Clifford algebra. 
Also M is a matrix. 

The usual Dirac equation where the matrices r y a fulfill a Clifford algebra and where 
M is proportional to the unit matrix, can be derived from the additional assumptions (v) 
uniquenes of the null cones, (vi) two helicity states only, and (vii) uniqueness of the mass 
shell, compare [|], ^) . From these demands we arrive at a Dirac equation in Riemann-Cartan 
space-time, where the coupling to torsion consists in the axial part only. All these principles 
arc operational since they can be proven directly by experiments. 

Some basic features of quantum theory mentioned above have been questioned previously 
and subsequently been tested or estimated from some existing data on atoms, for example. 
One of these features is the linearity of quantum theory which is basic in our understanding of 
all quantum phenomena. A generalized quantum field equation including a non-linear term 
has been introduced and discussed by e.g. Shimony and Weinberg (7), and references 
therein. Shimony himself proposed a neutron interferometry experiment which subsequently 
was performed by Shull et al || giving a strong restriction on the strength of a hypothetical 
nonlinear term in the Schrodingcr equation. Also spectral data of the hydrogen atom have 
been used for this purpose. Another features which has been discussed was the conservation 
of probability J9| . The second part of assumptions (v) to (vii) manifests itself in a breaking of 
local Lorentz- and local position invariance which can be tested by Hughes Drever, red-shift 
and atomic interferometry experiments |l(]] . 

In this paper we want to question the first of the four basic principles (i) to (iv) underlying 
our basic understanding of quantum theory, namely the unique deterministic evolution. Here, 
unique deterministic evolution means that if a quantum state tj}{x) is prepared at a time to, 
then the state is uniquely determined for times t > to- This implies that the evolution of 
quantum states is described by an evolution equation which contains a first time derivative 
only, that is, -^tp = Aip, where A is some operator. 

In order to test this principle, we propose a generalization of the usual Dirac equation 
by adding a second time derivative which violates this principle. This modified Dirac equa- 
tion is used to calculate modifications of the propagation of spin-i-particles as well as the 
corresponding modifications of the atomic spectrum (The hydrogen spectrum can also be 
used as justification of the four-dimensionality of space-time on atomic scales [0]). Up 
to now, all the experimental results are well explained using the standard theory, that is, 
the usual Dirac equation together with quantum corrections. If everything is well explained 
using the standard theory, then the modifications of the results due to the modifications 
in the Dirac equation, can be only smaller than the experimental error. Therefore, all the 
modifications can be restricted by comparing the calculated effects with the accuracy for 
the various measured effects. - However, future observation of neutrinos and high energy 
photons from gamma ray bursts (GRB) may be capable to distinguish between the various 
models of the Dirac equation. Nevertheless, one has to keep in mind, that on cosmological 
distances the parameters may depend on the position and thus the effect we are looking for 
may be cancelled during the propagation over long distances. Therefore, these observations 
have to be complemented by laboratory experiments like spectroscopy. 

2 The model: A modification of the Dirac equation 

The unique deterministic evolution implies that the evolution equation for the quantum field 
is of first order in the time derivative. This means especially, that the evolution equation 
is an equation without memory. In terms of a system of partial differential equations this 
means that this system should be of first order in time as it is the case for the Dirac equation 
ihdtip = —ihea ■ 'Vtp + firac 2 ^. If an evolution possesses a memory, the time derivative has to 
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be replaced by an operator, for example, an integral expression: Btp = —ihcot ■ Vi/> + /3mc 2 ip 
with, for example, Bip = J B{t,t')ip{t')dt' . In the case that the kernel B(t,t') of that kind 
of equation possesses certain properties (it should be analytic), then one can expand that 
kernel resulting in a system of partial differential equations with an infinite sum of terms 
containing time derivatives of arbitrary order, X)o° a i^t'0- 

Therefore, if the quantum field does not evolve uniquely deterministic or if quantum 
theory has a memory, then the resulting field equation in these cases contains arbitrary high 
orders of time derivatives. In a first approximation, this may be modeled by adding to a 
conventional quantum field equation like the Dirac equation or the Schrodinger equation, a 
term containing a second time derivative: 

d h 2 d 2 

ih—ib = —ihca ■ VV> + (3mc 2 ip — e —-^ tp . (1) 

dt mc 2, dt z 

In order to make e dimensionless, we introduced a factor 1/mc 2 in the term containing the 
last term. Here we assume that e is constant, i.e. does not depend on time or position. For 
e = 0, the above equation reduces to the usual Dirac equation. It is clear that the last term 
in (Q) violates Lorentz covariance. 

This modified Dirac equation can be derived from a Lagrangian: 

C = -ih (ip + d t ip - d t tp + ip) 

1 h 2 

+ -ihc (ip + a ■ VV> - • atp) - mc 2 ip + /3t/; + e -d t ip + d t ilj (2) 

2 mr 

= Uh^daii-da^i)) -mS^-e-^dt^dtiP. (3) 
This implies, in particular, that we have a conservation law 

= dtp + V • j (4) 



with the probability density 



and a current 



h 2 

p = — e (dtip + ip — ip + d t ip) (5) 



hctp + ctif; . (6) 



The coupling to the electromagnetic field can, as usual, be introduced by means of the 
minimal coupling procedure: 

C = —ih {dt + ie(f))ijj — (d t — iecj))i}} + if) 

h 2 

-mc 2 tjj + l3ip + e -{d t - ie(j>)tp + (d t + ie<j))xj; . (7) 

mcr 

The corresponding modified Dirac equation is 



ihdtip = -ihca • V - —A tp + f3mc 2 i> + ed>+ e ~{d t - iediYtb . (8) 

\ he J mc z 

3 Plane wave solutions and neutrino propagation 

It is not difficult to present an exact plane wave solution for ([!]). Inserting the ansatz 
ip = exp {i(Et — p ■ x)) a into (||) gives Ea = {a ■ pc + fimc 2 + ^sE 2 ) a, or[] 

7° (E - -^E 2 ) + 7 a p a c - mc 2 ) a . (9) 



1 Indices with a hat run from 1 to 3 and, for example, is represented by p. 
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The corresponding dispersion relation reads 

V mc z / 
which possesses the four solutions 



Eg 



4 2) 



i 2 c 4 =F 4emc 3 \/p 2 c 2 + m 2 c 2 
2e 

+ \J m 2 c 4 ± 4emc 3 -\/p 2 c 2 + m 2 c 2 
2e 



(11) 



For small e this reduces to 



2 

E { ± ] = ±c^p 2 +m 2 c 2 + e ' — 2 1 P 



mc 

TO 



_( 2 ) TOC 2 , ^5 5-5 ( 2 P 2 

= ± cyp + 7rr<r — e | toc - 1 - 



That means that for e — » the solutions E±^ reduce to the well-known solutions. The 
other two solutions are new and diverge for small e. However, this large quantity will drop 
out by considering energy differences. The ± in the solutions for the energy corresponds to 
positive/negative energies, the (1), (2) corresponds to the two solutions which come up with 
e ^ 0. 

In the high energy limit m, — > we get {p = \p\) 



(12) 



2 3 / 2 

„m toc / o V 

EX' = ±cp ± — — + e toc 2 + — 
zp \ m 

12) mc 2 , , m 2 c 3 ( n p 2 
EX = ± cp ± — e toc 2 + — 



(13) 



In the low energy limit p — > we get instead 

E± = ±mc ± h e toc H 

2?7l \ TO 

(2) TOC 2 2 p 2 / 2 p 2 

E± = ± toc ± e toc H 

e 2to \ to 

The exact expression for the group velocity reads 

(«£>)* = ± ^ 



(14) 



,/p 2 + m 2 c 2 J 1 =F 4-r-vV +?" 2 c 2 



(<4 2) r = ± 



(15) 



^/p 2 + m 2 c 2 J 1 ± 4^^p 2 + to 2 c 2 



Therefore, for positive energies, particles with the same momentum but with different direc- 
tions propagate with different velocities (the same is true for particles with negative energies). 
This property can be used for a comparison with data from neutrino propagation. 
For small e we get from dl5| ) 

i^f = ± /2 P&C 22 + 2£ - 

a/P 2 + to 2 c 2 m 

(v?f = ± , P& ° -2e^, 
Vp 2 + to 2 c 2 m 



for m — > 



, (l)xa i PS ^ ™ 2 C 3 Pa 



P 2p 3 "i 



/ (2)xa , ^ TOC p & 

p 2p J TO 
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and for p^ — > 

(„£))« = ± ?± + 2e El 

(18) 







+ 2e^ 


m 


m 


(«i 2) r = 


± Pa 

m 


m 



For a comparison with data from the propagation of neutrinos, which may be produced 
in connection with GRBs, we use (|l7|). We compare the arrival time of neutrinos with the 
arrival time of light over a distance of I = 10 10 ly. If the neutrinos and the photons are 
produced during the same event, and if we take the mass of the neutrinos to be 1 cV and the 
momentum p = 10 5 GeV, then we get as difference of the time-of-arrival 

M = ---^n-(^±2e-?-) « (4x 10- 9 ±6.4el0 31 ) sec. (19) 
c v ' c V °P mc J 

The first term can be neglected compared to the second one, so that we get |At| =Je|6.4 x 
10 31 sec. Asssuming a temporal structure of the source of about a millisecond and 
assuming a null-result, then we can get from observations of the propagation of neutrinos 
and of photons the estimate 

|e| < 1.6 x 1CT 35 . (20) 

Thus, neutrino observations in the future have the potentiality of high precision determination 
of the parameter e. Any |e| larger than that given by ( pp| ) should be detectable by this means. 

In quantum gravity theories, |e| is proportional to the ratio of the Planck length and 
some intermediate length, e = nl p /L where k is assumed to be of the order 1 If we take 
L = h/p, then, in terms of ft, the above estimate means \k\ < 1.6 x 10~ 21 which certainly 
is in contradiction to the assumption that k is of the order 1. From this we conclude, that, 
if the Dirac equation contains an additional quantum gravity induced term with the second 
time derivative of the neutrino field, then this term should be observable in the future by 
comparing neutrino propagation with photon propagation. 

However, from the derivation of the modifications of the Dirac equation [gj, the parameter 
e or k may be constant over the scale L of the "weave" states only. Therefore, it may be 
possible that the effect, as we calculated it, may not occur due to an averaging to zero over 
larger distances. Consequently, it is necessary also to perform tests which take place on a 
small scale only. One kind of such tests is atomic spectroscopy what we are going to discuss 
below. 

4 The non— relativist ic limit 

4.1 The non— relativistic field equation 

First we calculate the modified Pauli equation corresponding to eqn (Q). In order to do so, 
we first subtract the rest energy from the wave function by means of the substitution 

?A = e -^ 1+£ ) mc2 V (21) 
resulting in an elimination of the rest mass in one part of the wave function. This gives 

(l + 2e(l + e))iHd t ip' = -ihea ■ f V - tp' + (fi - l)mc 2 ip' + e<jyip' 

-e-^-idt + ie^'. (22) 



With the projection operators P± := ^(1 ± (3) we define the 'large' and 'small' parts of the 
wave function: V± = P±i>' ■ Multiplication of (|l|) with P + and P_ gives the two equations 

(1 + 2e(l + e))ihdt^+ = -itica • V - —A )ip'_ + e#+ - e ? (ft + ie<j>) 2 ^ + (23) 

\ nc J mc z 

(1 + 2e)ihd t tp'_ = -ihca. ■ (v - j^Aj + e< ^- ~ 2mc2 ^'- 

-e^idt+ie^L- (24) 

As usual, we assume that in the second equation the energies ihdtip'_ and afiip- are small 
compared with the rest mass term mc 2 ip'-- Therefore, we approximate 

*'=-2^ ifa '( v -i A K m 

Inserting this into the first equation ([23]) gives 





h 2 / 

2m V 


V 






mc 
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(ft + ie^) V+ , (26) 

I ILL, //tU- 

with S = I ^ ° ) and B = V x A. Here <r arc the three Pauli matrices. After division 



rr 

through 1 + 2e(l + e) we get 
*2 / 



n 2 / ?> \ p /> 2 

W ; = -— (V - -A) tf + + eW + + ^ " B - ^(ft + ie *)V + , (27) 

where we absorbed the e-factor in a redefinition of mass and charge: m* := (1 + 2e(l + e))m, 
e* = e/(l + 2e(l + e)). 

As our result we get a modified Schrodinger equation for a two-spinor i/j 

8 h 2 ( ie \ 2 e fr 2 

lh 7T^ = — o — * V " "T A ^ + + —<t B- e—^^dt + ie<f>) 2 iP . (28) 
at 2m \ nc / mc m c 

This equation can also be derived from a Lagrangian 

C = l -h{il}*{d t ~ie(t>)i}-{d t +ie(t))^*i:) (29) 



2m* V fic 



a) V* • fv - V-' - (ft - ie0)^*(ft + ie$i& . 

/ \ nc J m cr 



Equation ( pq ) also possesses plane wave solutions whose energies are given by (|l4]). We 
also have a conservation law 

|p+Vi = (30) 



with the probability 



and the current 

2m* 

Like in the Klein-Gordon equation there seems to exist the possibility to get negative prob- 
ablitics. However, using the Schrodinger equation, we get for the probability 

1/9 d 

p = ip*ip - e— j I -ih—<ip*ip + %p*ih—ip 

= ^ + c _L J C^-A^ + V*^A^y (33) 
m*c z \2m* 2m* / 

This quantity is strictly positive if e < 0, and for e > this is positive if |e| is small enough. 
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4.2 Modifications of the energy levels 

It is possible to calculate the energy levels of the hydrogen atom exactly. In order to do so 
we choose A = and cf> = <p(r) as the usual spherically symmetric electrostatic potential: 

d H 2 1 ( \ 2 



dt 2m* m*c 2 \ dt 

We asume a stationary solution, then ihd t ip = Eip, and we get from the above Hamiltonian 

(e - e-^E 2 ) $ = -^Ai> + (e* - 2e^e) # + e-^e 2 ^ 2 ^ . (35) 
\ m*c z J Ira* \ m*c z / m*c z 

With A = ( r2 -§^) — piL , where L is the angular momentum operator with the eigen- 

value equation L Y™ = 1(1 + 1)Y™ with I = 1, 2, 3, . . . . With a splitting of the wave function 
into a radial and an angular part ip = i?(r)Yj m (■# , ip) we get the radial part of the wave 
equation 

Jl2 ; 



2m* \r 2 dr \ d 



e*-2e kE 4>R + e T e 2 ffi. (36) 

m*c / m*c 2 

With the explicit expression for the electrostatic potential (f> = —e/r we get by multiplication 
with 2m* /h 2 



d 1 „ 2d 2?n* / 1 



= —R+-^-R+-—[E-e—-rE 2 )R 



dr 2 r dr h 2 \ m*c 2 

2m* f „ „ e „\ e „ 1 /,„ 2e 



V m c / r r z \ nrc J 

Since the r-dependence is the same as for the usual hydrogen atom, this equation can be 
solved using the standard scheme: With 

where a = e 2 /he is the fine structure constant, we get for the radial part of the wave function 

^ + 2 -± + E + 1-^)r = 0. (39) 
dr £ r dr r r z J 

We introduce dimensionless coordinates by 

r' = -, -S = ^j, q' = qro (40) 
r 4rg 

and get with R'{r') := R(r(r')) 

d 2 2d 1 q' 1 



, )R' = 0. (41) 
dr' 2 r' dr' 4 r' r ' 2 J 

We also introduce a new variable f(r') through 

R = e -i r 'r'~<f(r') (42) 

with a parameter 7 which will be specified later. We get the following equation for the 
function /: 

..^.^.^(iijh!^.,.^,,. (43 ) 
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In order to solve this equation we specify the value of 7 by the requirement that the term 
~ 1/r' should vanish: 

7(7 + 1) -£ = 0. (44) 

This gives the two possibilities 



7 ± = ~\ ± >fi^\ = ~\ ± y (l + 5) 2 + 2«* 2 , (45) 

and from (^3|) the differential equations 

z/" + ( t 9 ± -z)/'-/3±/ = 0, (46) 

with 



e + - (47) 



/?± = 1± + l- q > = l±j£+l- q >. (48) 



Eqn. (|46|) is the confluent hypcrgcomctric differential equation with the solution |1 



!/=0 ' 

which is appropriate for our problem. 

It is clear that, in order to get no infinite terms, $± is not allowed to be a negative 
integer: $± ^ — f , —2, . . . , which is fulfilled if e 7^ and |e| < 1. For e = 0we cannot use the 
solution In addition, if the sum in (49) does not terminate, then the solution diverges 
for large r faster than exp(ir') which leads to non-normalizable solutions. The condition 
for a termination of the sum is /3± G Z~ , or 

&± = -k, k e N . (50) 

With (f48|), (Eoh, and (p8h we get four energy eigenvalues E 



(1) m*c 2 K 2 ± ~2eaa*^ m*c 2 K ±\J K ± + 2eaa*{l - 2jj 

E± ~ ~2k ~Rg - 2ea 2 + ~e Kg - 2ea 2 (51) 



(2) m*c 2 K% -2eaa*^ m*c 2 K ±\f K ± + 2eaa*(l - 2^) 

S± ~ ~2e Kg - 2ea 2 Te Kg - 2ea 2 (52) 



with K± = k - i =p y + l) 2 + 2«* 2 and a* = e* 2 /he. To first order in e we get 
m) „ 2 / 1 a 2 1 4 5 - 8fc + 21 

E + = m c [I + T (I-IT7F - a £ 4(i - k + o 4 d + 20 + 



- - V U + T(fcT7F + - 4 4 (fc + 04 (i + 20 + 

< = ~ m c ( t (i - fc + if - ea 4(i - 3b + mi + 2/) + 0(6 > ) (55) 



,(2) ™*„2^ 2 1 , _ 4 3-8fe-2Z 



E - = ~ m * c U(feT7F + £a 4 (fc + oni + 2O +0(£ ) J- (56) 

The solution is then given by 

fl(r) =i2(r(r / )) = -R'(r') = Ae~? r ' r'i f(r') = Ae~i r ' r n± f(-k,d ± ,r') (57) 
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where A is a normalization constant and where all parameters depend on I. 

For e ^ all these energy values are well defined. Even for very small e the first two 
energies ( |53"| , j54| ) are valuable solutions, too, since only energy differences are observable and 
thus the first term drops out. However, there are two reasons which justify to drop the first 
two solutions: (i) Except the first term mc 2 /e, the upper two sets of energy levels ( [53| , [54| ) 
are the same as the lower two sets ([35 56 ) up to sign. For the first two sets of energy levels 



the continuum is below the discrete spectrum. By postulating that all particles fall into the 
lowest energy level, then all atoms will fall into the continuous part of the spectrum which 
has never been observed. (One also can postulate that all particles want to go to the highest 
energy level. But due to the symmetry of the two sets of spectra, this will give the same 
answer. Therefore we do not consider the upper two sets by convention.) (ii) We want to 
describe small modifications of the known energy levels given for e = 0. Owing to these 
reasons, we keep the last two sets of energy levels. 

Defining for E^ the principal quantum number n := 1 — k + 1, then the energy levels are 
given by 



1 . ,/J3 

2n 4 n 3 (l + 2Z) 



4 2) = -™'*=r liW A- -^TT ) + 0(e*) ) (58) 



and if we define for E^ the principal quantum number n := k + 1, then the energy levels are 
h(2) _ m , j" 2 ( 1 , 2 ( 3 



where mc 2 a 2 /2 is the Rydberg constant Ry. It is remarkable, that in both cases we get the 
same scheme of energy differences. Therefore, at this order we cannot distinguish between 
energies E+ and E-. 

There are two differences between this result and the usual spectrum of the hydrogen 
atom: (i) No degeneracy with respect to I. (ii) Additional 1/n 4 - and l/n 3 -tcrms. From ( [58| ) 
we get for the Lyman-series with n = 1 and I = as ground state, for example, the energy 
differences 

AE — - Ry { 1 -^ + {li 1 -^)-^ 1 - )) + °&) • ( 6 °) 

The additional terms 1/n 4 and n 3 (l + 21) modify the structure of this series. The ionization 
energy is given by n — > oo and is 

^ionization = Ry U ~ \<iO? + 0{e 2 )^j . (61) 

From numbers related to AE'Lyman or -^ionization one can draw estimates on the value of e. 

Since the accuracy SAE/AE for recent measurements, see for example |l4|], is of the order 
of 10~ 13 which agrees completely with the conventional theory, we can conclude that the 
relative deviation from the usually calculated energy has to be smaller than this uncertainty: 



SAE 
AE 

This means that 



5 3 4 1 

= ea 2 \~ 2 ~^ 1 + ^ < 1(T 13 . (62) 



1-A 1 



* < iO" 13 ^^ T^-\ = lO- 13 -^-^ « 7 x 1(T 10 . (63) 



5 3_ j 4_ I n 2 a , 3 

' 2 2ri 4 "•" «3 I " ° "T 8 



In terms of a k as introduced after Eqn. ( |20| ) we have the estimate k < 7 x 10 4 which is outside 
the assumption that k is of the order unity. Therefore, the accuracy of atomic spectroscopy 
is still at least five orders of magnitude too small in order to be able to detect any influence 
of quantum gravity on atomic levels. 
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5 Conclusion 



We have shown that an additional term in the Dirac equation containing a second time 
derivative, as it is motivated from quantum gravity, influences neutrino propagation and 
atomic spectroscopy. While neutrino propagation may be capable to 'see' this additional 
term, the accuracy of spectroscopy has still to be improved by five orders of magnitude in 
order to be sensitive to this term. 

A clear difficulty with the hydrogen spectrum is that, precisely, it cannot be calculated 
exactly with an arbitrary accuracy with the present state-of-the-art. There is a number 
of corrections to the Dirac solutions: recoil, QED, finite nucleus size, sec |l6| for a review. 
These corrections scale with 1/n 3 and amount to a few kilohcrtz. They are very difficult to 
calculate and there are still some discrepancies between theoretical results. Fortunately, some 
combinations of frequencies are independent of these corrections at their leading order and can 
be used for higher accuracies. Concerning the measurements themselves, there is presently 
a very rapid evolution towards much higher accuracies. For example, it was shown recently 
that, because of the extreme regularity of the frequency comb generated by femtosecond 
lasers over a very wide spectrum [[j~7| , they could be used to compare frequencies of oscillators 
which differ by several orders of magnitude. This provides a way to compare many transition 
frequencies of the hydrogen atom between themselves and with microwave clocks, with the 
potential accuracy of the cesium fountain clock which is presently 10~ 15 and should improve 
quickly by another order of magnitude. Also, the techniques to interrogate narrow transitions 
of cold atom hydrogen by sub-Doppler methods or atom interferometry have improved very 
significantly either in cold thermal beams (Hansch and coworkers in Garching, Biraben and 
coworkers in Paris) or in clouds generated from Bose-Einstcin condensates (Kleppner and 
coworkers at MIT). A subkilohcrtz linewidth is presently achieved for the 1S-2S two-photon 
transition [ p~7[ and could still be reduced by one or two orders of magnitude in the near 
future. The hydrogen atom is thus potentially a universal clock by itself covering the full 
spectrum from UV to microwaves (hydrogen maser). 
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